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ABSTRACT 

The radio and X-ray observations of SN 1993J during the first year can be 
consistently explained as a result of interaction of the expanding ejecta with a 
circumstellar medium. The density of the circumstellar gas can be deduced from 
the free-free absorption of the radio emission and from the X-ray luminosity. During 
the first two weeks, both sets of observations indicate a mass loss rate of ~ 4 X 
10 - ^ Mq yr~l for a wind velocity of 10 km s - ^. The subsequent radio and X-ray 
observations indicate a density gradient oc r~ s , with 1.5 ^ s ^ 1.7, as opposed to the 
r gradient expected for a steady, spherically symmetric wind. This may either be 
caused by a variation of the mass loss rate from the progenitor system, or by a non- 
spherically symmetric geometry. To explain the properties of the X-ray emission, a 
steep density gradient in the ejecta is needed. During the first months most of the 
observed X-ray emission originates from the circumstellar shock, which is adiabatic, 
while the reverse shock is radiative. To avoid excessive Comptonization in the X- 
ray range collisionless heating must be ineffective. The soft X-rays observed at 220 
days probably originate from the reverse shock. The ionization and temperature 
structures of the circumstellar gas are calculated; we find that the temperature 
is in excess of 10^ K and the medium is nearly completely ionized by the shock 
radiation after the formation of the shocks. Preacceleration of the circumstellar gas 
by the radiation from the outbreak can explain the observed high velocity for the 
circumstellar N V and Ha lines. The high luminosity of the lines indicate that the 
circumstellar medium close to the supernova progenitor had a complex structure. 

Subject headings: circumstellar matter - 7-rays: general - radio: general - 
supernovae: individual (SN 1993J) - X-rays: general 

1. INTRODUCTION 

SN 1993J was discovered on March 28.86 1993 (Ripero & Garcia 1993). The 
first detection was on March 28.30 at magnitude 13.4, while a photograph taken on 
March 25.6 showed no supernova brighter than 17 magnitude (van Driel 1993). The 
explosion date is thus confined to between March 25.6 and March 28.3. Based on the 
apparent rapid rise around March 28.0 and the roughly 0.2 day for the shock wave 
to traverse the star (T. Shigeyama, private communication), we adopt March 27.8 
as the time of the explosion. There is evidence for the progenitor object (Perelmuter 
1993), but its exact nature is still uncertain. Aldering, Humphreys & Richmond 
(1994) have suggested that the photometry is consistent with two supergiants of 
spectral classes early B and early K. Based on the early light curve, Nomoto et 
al. (1993), Podsiadlowski et al. (1993), Ray et al. (1993), Bartunov et al. (1994), 
Utrobin (1994) and Woosley et al. (1994) interpret the progenitor as a star of initial 
mass ~ 15 Mq, which has lost all but 0.1 — 0.9 Mq of its hydrogen envelope. They 
propose that Roche lobe overflow in a binary system is responsible for the loss of 
the hydrogen envelope. A substantial fraction of this mass may have been lost from 
the system. An alternative scenario based on a single massive star with an initial 
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mass of 30 M©, which has lost most of its hydrogen envelope in the form of a wind, 
has been proposed by Hoflich, Langer & Duschinger (1993). In both scenarios one 
expects a substantial circumstellar density. 

Direct evidence for the presence of circumstellar gas comes from observations 
with IUE and optical telescopes. In particular, IUE observed a strong N V line at 
1238 - 1242 A (Sonneborn et al. 1993; Fransson & Sonneborn 1994) with a width 
extending up to ~ 1000 km s - ^. From the first observations a few days after the 
explosion, the flux of the line decreased by a factor ~ 100 on a time scale of about 
a week. The width implies that it must arise outside the ejecta. The high degree 
of ionization and fast decay constrain it to be formed close to the supernova. This 
is in contrast to SN 1987A where the emission increased on a time scale of ~ 400 
days as a result of a light echo in the ring (Lundqvist & Fransson 1991). In the 
optical, narrow lines of Ha, He II A 4686, and possibly [Fe X] A 6374, [Fe XI] 
A 7892 and [Fe XIV] A 5301 were observed (Filippenko & Matheson 1993a; Benetti 
et al. 1994; Cumming et al. 1994). The Ha line was resolved with a velocity similar 
to that of the N V line (Cumming et al. 1994). The Fe lines had a considerably 
smaller FWHM of only ~ 50 km s — On April 13, the narrow lines were no longer 
observable (Filippenko & Matheson 1993b; Cumming et al. 1994). 

X-ray emission was first detected by ROSAT on April 3.4 in the 0.1 - 2.4 keV 
band (Zimmermann et al. 1994) and on April 5.25 by ASCA in the 1-10 keV band 
(Tanaka et al. 1993). OSSE on the Compton Observatory detected the supernova at 
^ 100 keV (Leising et al. 1994). The early appearance suggests that the emission 
has a circumstellar origin, in contrast to the Comptonized emission from the down 
scattered 7-rays from the radioactive decay. 

Radio emission has been observed with a first possible detection by the VLA 
at 1.3 cm on April 2.30. Radio turn-ons have since been observed at subsequently 
longer wavelengths, consistent with free-free absorption in circumstellar gas (Weiler 
et al. 1993). On April 5.7 it was clearly seen with the Ryle telescope at 2 cm with 
a flux of 3 mJy; no emission was seen the previous day (Pooley & Green 1993). 
Subsequent observations at this wavelength showed a roughly linear rise up to day 
60. The supernova has also been observed at 3 mm (Phillips & Kulkarni 1993a, b). 

The presence of radio and X-ray emission together with narrow emission lines is 
not new for Type II supernovae. For example, both SN 1979C and SN 1980K showed 
strong radio emission, turning on some time after the optical outburst (Weiler et 
al. 1986). In addition, SN 1980K was detected as an X-ray source by the Einstein 
Observatory at a luminosity of ~ 2 X 10^ erg s - ^ (Canizares, Kriss, & Feigelson 
1982). All these features have been interpreted as a result of interaction of the 
supernova ejecta with its circumstellar medium (see Chevalier 1990 for a review). 
Both these and several other Type II supernovae have been detected several years 
after the explosion in the optical, long after radioactivity has decayed away, which 
is consistent with circumstellar interaction (Chevalier & Fransson 1994; hereafter 
CF94). In this paper we show that the observations of SN 1993J fit into this picture, 
although there are some deviations from the 'standard' model. Throughout this 
paper, we take a distance to M81 of 3.6 Mpc (Freedman et al. 1994). 
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2. IMPLICATIONS OF THE RADIO LIGHT CURVES 



Assuming that the progenitor was a late type supergiant, it is expected to have 
had a slow moving dense stellar wind. If the mass loss rate is M, the wind velocity 
is v w and the mass loss is spherically symmetric, the number density at radius r 
is n cs = M I (Att pmjjv w r^), where p is the mean atomic weight in amu and mjj is 
the mass of a hydrogen atom. For most of our discussion the density is the relevant 
quantity, and only the ratio M /v w enters as a parameter. For typical red supergiants 
the wind velocity is 10 — 20 km s - ^ and the mass loss rate is 10 - ^ — 10 - ^ Mq yr - ^. 

The observed radio turn-ons of previous Type II supernovae have been well 
described by a model where the turn-on is a result of the decreasing free-free 
absorption of the ionized stellar wind as the supernova expands (Chevalier 1982b). 
Although neutral before the outburst, the radiation from the photosphere at the 
time of shock break-out and from the shock region ionizes and heats the wind to a 
high temperature (Lundqvist & Fransson 1988). The interaction of the ejecta and 
stellar wind results in a shock wave propagating into the stellar wind, and a reverse 
shock into the ejecta. 

Hydrodynamical calculations and analytical arguments show that the ejecta 
density is well described by a power law in radius (or velocity). A similarity solution 
for the shock radius can then be found (Chevalier 1982a, b). We write the ejecta 
density as a function of velocity, v, as p sn = po(t/tQ)~^(v /vq)~ u , where pq is the 
density at velocity vq and time £q- For compact progenitors, the power law index 
is expected to be in the range 7-12 (e.g., Chevalier & Soker 1989; CF94). An 
extended red supergiant may, however, have a considerably steeper gradient. The 
only two independent quantities are then pq^Vq and M/v w . If both shocks are 
radiative, as may be the case shortly after the explosion, a thin shell approximation 
is applicable. Below we show that a circumstellar density variation different from 
p w oc r is indicated from the observations. We therefore write 

M fr \ s , s 

Aitv w r% \ r J 

where r is a reference radius, corresponding to the mass loss rate M/v w . In the 
following we use r = 10^ cm. The maximum velocity of the ejecta, V = R s /t, is 
then given by 



V 



4tt (3 - s)(4 - s) P Q 4 vfi v w rl~ s 
(n - 3)(n -4) M 



l/(n-s) 

t -(3-s)/(n-s)_ ( 2 _ 2 ) 



More general cases differ from this approximation by less than ~ 20%, and are 
discussed in CF94. The outgoing circumstellar shock has a velocity v cs = £ V(n — 
3)/(n — s) and the reverse shock v rev = (3 — s)V / (n — s), where £ is the ratio of the 
forward shock radius to the reverse shock radius. Because V is in principle directly 
observable as the maximum ejecta velocity, we will use it as a parameter in this 
paper, rather than pq. 
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The free-free absorption optical depth of the unshocked ionized gas at wave- 
length A is given by 

T//K0 , 6 ^)Vvr-(^p (,3, 

where M_5 is the mass loss rate in units of 10 - ^ Mq yr~\ v w \ the wind velocity 

in units of 10 km s~\ T5 the wind temperature in 10^ K, and V4 the ejecta velocity 

at the reverse shock in 10^ km s — We have evaluated the Gaunt factor for A = 2 

cm, a number density ratio n(He)/n(H) = 0.10, and T = 10^ K, and assuming 

£ = 1.3. The data of Van Dyk et al. (1993) on 22.5 April indicate that r = 0.5 

at 2 cm at that time. The shock velocity is not directly measurable from the 

spectra, but is expected to be somewhat larger than the highest ejecta velocity. 

The blue edge of the Ha line showed a velocity of ~ 19000 km s - ^ on 1993 

April 11 (e.g., Lewis et al. 1994). Earlier spectra show evidence for even higher 

velocities. The shock front is more directly observable at radio wavelengths using 

VLBI techniques. The source was resolved by Marcaide et al. (1994) on 1 May, 

yielding a velocity of (2.4±0.15) X 10^ km s - ^ for a uniform disk model, and by Bartel 

et al. (1994) on 27 April, 17 May, and 27 June, 1993, yielding velocities in the range 

18000 - 19000 km s" 1 with 5-10% uncertainties. We will use V = 20000 km s" 1 

on day 10 for most of the paper. When this is substituted in equation (2.3), and 

taking limb darkening into account, which reduces the mass loss estimate by ~ 20%, 
/ 3/4 

one obtains M_5/u u ,i 3.8 Tg . The temperature of the unshocked wind close 
to the shock depends on the properties of both the stellar wind and the expansion 
rate of the supernova, but is typically ~ 10^ K (see § 4.1). We therefore find that 
M_$/v w i 4, within approximately a factor of 2. This is an order of magnitude 
larger than that determined by Panagia, Van Dyk & Weiler (1993). The reasons 
for this are probably that Panagia et al. use a wind temperature of 10^ K and an 
expansion velocity of 10^ km s - ^. 

Using the model by Chevalier (1982b) and including limb darkening effects by 
the wind, we have modeled the radio light curve at 2 cm obtained by Pooley & 
Green (1993). We assume that the magnetic field and relativistic particle energy 
densities scale as the shock ram pressure, pv^ . In figure 1, we show our best fit for 
s = 2, n = 30, spectral index 7 = 0.7 (flux oc v~^) and r(2 cm) = 0.5 on 22 April, 
1993. The fit is considerably worse than for earlier, "normal" Type II supernovae 
(e.g., Chevalier 1984; Weiler et al. 1986). The differences between the model and 
the observations are especially large for the rising part of the light curve. We have 
therefore explored varying some of the parameters with respect to the standard 
model. 

The rising part of the light curve is determined mainly by the exponential de- 
crease of the optical depth with time, F\ oc e~ Tx ^ r ^\ and is not appreciably affected 
by assumptions about the energy density of the magnetic field and relativistic par- 
ticles. From equation (2.3) and with T e constant (§4.1), the linear shape of 
the light curve therefore indicates a density variation different from the standard 
p w oc r~ 2 . 
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The optically thin part of the light curve is given by 

F oAK A^r{ 3+1 - [6+1 -* (3+1 »^}. (2.4) 

Type II supernovae typically have 7 0.8 (e.g., Weiler et al. 1994). Therefore, a 
nearly constant flux may be produced by a radial wind with p w oc r - ^. We find 
that good fits are indeed obtained for p w oc r - ^. This is shown in figure 1 as the 
dashed line for n = 30, 7 = 0.7, n(He)/n(H) = 0.1, T5 = 1, and M-^/v w \ = 1.35 at 
10^ cm. The parameters were chosen to fit observations also at other wavelengths. 
In figure 2 we show observations by Phillips & Kulkarni (1993a) and Weiler et al. 
(1994) together with the p w oc model. 

The epochs of the radio turn-ons at different wavelengths are sensitive to the 
value of n. A lower n than n = 30 would cause a larger spread in the turn-on 
times at the different wavelengths than in figure 2. Large values of n were also 
obtained for SN 1979C and SN 1980K (Chevalier 1984; Weiler et al. 1986). The 
radio evolution in the standard model depends on the assumption that both the 
particle density and magnetic energy density scale as the ram pressure. If, e.g., the 

0. 

magnetic field decreases slower than B oc pv , this affects the turn-on in roughly 
the same way as a large value of n. The departure needed from a simple scaling is 
a factor of ~ 2 in the for the p w oc model over a factor of ten in radius to 

reduce the value of n in a best fit model from n = 30 to n = 15, though the fits are 
not as good as in figures 1 and 2. For lower values of n the fits are not compatible 
to the observations. A change in the relativistic particle density has roughly the 
same effect on the radio light curves as a change in the magnetic field. 

The observations at 3 mm (Phillips & Kulkarni 1993a, b) show a drop in lumi- 
nosity after ~ 100 — 150 days which is consistent with the law extending only 
out to ~ (2 — 3) X 10^ cm. At larger radii the observations are consistent with an 

law. In §5, we discuss further the implications of the fact that an wind 
fits better to the radio data than steady, spherically symmetric mass loss during 
the first ~ 100 days 

3. X-RAY EMISSION 

X-ray emission was first detected by ROSAT on April 3.4 (day 7) with a 

luminosity of 2.94 X 10^ erg s - ^ in the 0.1 - 2.4 keV band (Zimmermann et al. 

1994). Subsequent observations showed that the luminosity decreased by ~ 44% 

from April 3.4 to May 4. It was also detected on April 5.25 by ASCA in the 1 

- 10 keV band with a luminosity of 5 X 10^ erg s - ^ (Tanaka et al. 1993). For 

a thermal bremsstrahlung spectrum the ASCA team found kT > 10 keV and a 

column density ~ 10 cm - . The latter is somewhat larger than deduced from the 

on 

ROSAT observations, ~ 5.9 X 10 cm , which were consistent with the galactic 

on 

absorption in the M 81 direction, ~ 4.3 X 10 cm" . 

Observations using OSSE on the Compton Observatory have revealed SN 1993J 
as a source of hard photons between 50 - 150 keV (Leising et al. 1994). During the 
first two observations on days 9-15 and days 23 - 36 the intensities at 100 keV were 
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1.82 x 10~ 3 photons cm -2 s _1 MeV -1 and 0.89 x 10~ 3 photons cm -2 s _1 MeV -1 , 
respectively. A third longer observation on days 93 - 121 gave only an upper limit. 
With photon spectral indices, -2.3 and -2.2, respectively, the luminosities between 
50 - 150 keV are 5.5 X 10^ erg s _1 on day 12 and 3.0 X 10^ erg s _1 on day 28.5. 

The X-ray emission may originate either in the reverse shock or the circumstellar 
shock. In the thin shell approximation the temperature of the wind shock is 

T cs = 2.27 x 10 9 (i 8 ^ ~ 3 ^ Vl K, (3.1) 

[n — sy 

while the reverse shock temperature is 

(3-3) 2 
(n-3) 



T r ev — ( _ Test (3-2) 



assuming electron-ion equipartition in both cases. The mean mass per particle 
(electrons and ions) in a.m.u. is p s = [1 + 4 n(He)/n(H)]/[2 + 3 n(He)/n(H)]; for 
solar abundances p s = 0.61, while for n(He)/n(H) = 1, as may be the case for SN 
1993J (Baron et al. 1994), p s = 1. If electrons and ions are not in equipartition the 
ion temperature is a factor ~ 2 higher than the above estimates. 

The density behind the circumstellar shock is p cs = 4 p W} where p w is given by 
equation (2.1). The density ratio of the reverse and circumstellar shocks is given by 

p ~ (n - 3)(n ~ 4) p (3 3) 

The swept up mass behind the circumstellar shock is 

and that behind the reverse shock 

M rev = ^4 ( 3 - 5 ) 
(4 - s) 

For n = 6 — 30, the temperature behind the wind shock is <; 10 9 K, while for the 
reverse shock the temperature is 10^ — 5 X 10^ K. The time scale for equipartition 
between electrons and ions is 

* e -,«2.5xlO«( I |^) 15 ( I ^)" 1 S (3.C) 

(Spitzer 1978; Stepney 1983). At 10 days n e < 3 x 10 8 cm" 3 and T e < 3 x 10 8 K 
behind the reverse shock, so t e _{ ^ 8.5 days; the reverse shock is marginally in 
equipartition, even without collisionless heating. The ion temperature behind the 
wind shock is (1 — 2) X 10^ K for V\ = 2, and the density is a factor <; 4 lower than 
that behind the reverse shock. Electron-ion collisions are therefore highly ineffective, 
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and the electron temperature may be much lower than the ion temperature, unless 
we invoke efficient plasma instabilities heating the electrons collisionlessly in the 
shock front. 

The higher density behind the reverse shock means that the X-ray emission 
below 10 keV is dominated by the reverse shock (see below). A possible exception 
may occur if inverse Compton scattering of photospheric photons by the relativistic 
electrons in the shocked region is important. However, the observed flat spectrum 
(Tanaka et al. 1993; Zimmermann et al. 1994) is incompatible with a spectrum dom- 
inated by Comptonization. Below we use this to constrain the electron temperature 
of the wind shock. 

For the reverse shock to dominate the observed flux, absorption by the postshock 
gas cannot be large, as is inferred from the ROSAT observations during the first 
~ 100 days. A cool, absorbing shell will be present if the reverse shock is radiative. 
Therefore, for the observed X-ray flux to come from the reverse shock a necessary 
condition is that this shock is adiabatic. At temperatures above ~ 2 X 10 7 K, 
bremsstrahlung dominates the cooling, and the cooling time behind the reverse 
shock is t coo i = 1.7 x 10^ T^ e y/n rev s, or 

• -1 / ^ \ 4-s 

t coo i = i.6 x io° n ,_"- v ' ow 7 A , (^) v£-° 1 \ 




11.57 days 



In Table 1 we give C n s f° r different values of n and s, together with t coo { for 
the specific case M-^/v w \ = 2 and ^(lO days) = 2. It is seen that t coo \ de- 
creases rapidly as the density gradient gets steeper. According to equation (2.2) 
t cool oc M-(15+s 2 -9s+n)/(n-s) t (4n-ns-15+4s)/(n-s) ^ and cooling is therefore most 
important at early times. If cooling occurs early on, cool gas is likely to be present 
between the reverse shock and the contact discontinuity later, giving rise to an X- 
ray absorption much higher than observed (e.g., CF94). The importance of cooling 
therefore depends on the time of shock formation. To be conservative we evaluate 
t coo i one day after the explosion in Table 1. For M_§/v w \ ^ 2 and ^(lO days) ^ 2 
the reverse shock becomes radiative for n ^ 9. An earlier date would decrease this 
limit further. Equation (3.7) assumes that only bremsstrahlung is important for the 
cooling. In §3.2 we discuss the case when T e < 2 x 10 7 K, as is the case for n ^ 15 
(eqn. [3.14]). 

Summarizing, there are therefore two possible scenarios for the X-ray emission. 
Either both shocks are adiabatic, and most flux below ~ 10 keV comes from the 
reverse shock, or alternatively, the flux from the radiative reverse shock is absorbed 
by cool gas, and all the observed flux originates from the circumstellar shock. The 
latter case occurs for a steep ejecta density gradient. We discuss each of these cases 
separately below. 

In either case we can estimate the free-free luminosity from the two shocks from 
Ljg = Att f j jj{T t ){pjmii)^r^dr jjj(T^,) p) s /m'jj, where k denotes either the 
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reverse or circumstellar shock. The emissivity per unit mass is j ffPl m \ with 

jff(E) = 1.64 x 10" 20 ( 9ff T~°- b e ~ E / kTe erg cm 3 s" 1 keV" 1 , (3.8) 

where gjj is the free-free Gaunt factor, which includes relativistic effects, and 
C = [1 + 2 n(He)/n(H)]/[l + 4 n(He)/n(H)]. For n(He)/n(H) = 0.1, ( = 0.86, and 
for n(He)/n(H) = 1, ( = 0.60. Because of occultation by the ejecta, only half of the 
above luminosity escapes outward. 

As long as the emission in the ASCA and ROSAT energy ranges comes from 
the same region the ratio of the ASCA/ROSAT flux is given by 

LASCA _ A 1Q keV V 9ff(E) e~ E / kT c dE _ ,-7.9xl0 7 /T fl ( o ^ 



" Jo.fkeV 9ff(E) e-E/kT e dE 



within ~ 5% for 10^ ^ T e ^ 10^ K. The observed ratio is ^ 1.5, which sets a lower 
limit to the electron temperature of ~ 2 X 10^ K for the X-ray emitting gas. 

The density behind the reverse shock is large and the cooling time may become 
short before the contact discontinuity is reached. Lack of electron-ion equipartition 
in the circumstellar shock leads to an electron temperature much lower than the 
ion temperature. Cooling and deviations from equipartition modifies both the 
density, temperature and the emission from the shock, compared to the adiabatic 
similarity solution. We have therefore calculated the hydrodynamical evolution 
of the interaction numerically, using an explicit, Lagrangian code with artificial 
viscosity. Radiative cooling is included according to Raymond, Cox, & Smith 
(1976), as well as Compton cooling and heating by both the radiation from the 
shocked gas and the diluted photospheric flux. Above ~ 2 X 10^ K cooling of 
the reverse shock is dominated by free-free emission. We calculate this taking 
relativistic effects into account, including both electron-ion and electron-electron 
bremsstrahlung (e.g., Stepney & Guilbert 1983). The case that the electrons are 
assumed to be heated by Coulomb collisions only (Spitzer 1978; Stepney 1983) is 
referred to as T e = Tfj ou i. To examine the importance of collisionless heating we 
have assumed that T e = T{ on in some models. 

3.1. The adiabatic case 

We first discuss the adiabatic case. This may not be relevant for SN 1993J, but 
is of interest for other supernovae with shallow density gradients of their ejecta. In 
this case the OSSE flux is dominated by the circumstellar shock, unless the ejecta 
density gradient is extremely flat, while the reverse shock is responsible for the low 
energy flux. 

The lower limit of 10 keV for the reverse shock temperature at 10 days (Tanaka 
et al. 1993) implies in this case that n ^ [4.4 V\ (3 — s) + s]. For V\ < 2 

we find n ^ 15. Any departure from electron-ion equipartition decreases this limit 
further. As we found above, an adiabatic shock requires n ^ 9, so this constraint is 
consistent. 
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At ~ 1 keV, g ff « 1.87 Tg - 264 , where T 8 = (T e /10 8 K). The ROSAT luminosity 
from the reverse shock is therefore 

L rev (l keV) « 1.74 x 10 37 ( ^ ~ ^ ~ t V" 24 e" - 116 ^ 



3-2s ( 3 - 10 ) 
M_R V K 3 " 2s ( ..J, 1 erg s" 1 keV" 1 



4 I 11.57 days 



For a flat spectrum, the ROSAT observation at 7 days corresponds to 1.3 X 
"Lq39 er g g — 1 jj e y— 1_ With Tg ^ 1, as indicated by the observations, and V\ = 2, 
and assuming that the reverse shock is responsible for the X-ray emission, we find 
M_$/v w i = 2.9, 1.7, and 1.1 for n = 6, 7, and 8, respectively. With s = 1.5 the cor- 
responding mass loss rates measured at 10^ cm are M_^/v w i = 5.0, 2.9, and 1.9, 
respectively. The non-uniform density behind the reverse shock tends to decrease 
the mass loss rate by up to a factor two (see below). 

For T rev ^ 5 X 10^ K the reverse shock contributes only marginally to the OSSE 
luminosity. With flf//(100 keV) 0.32 (1 + 3 Tg) in the range 0.1 < Tg < 10, the 
luminosity at 100 keV is 



T cs (100 keV) « ( T e / 5 (1 + 3 T e9 ) e" 1 - 16 ^ f ^ 

(3 - s) \ v wi 



T/ 3-2s ( fj days \ 3 2S -1 , ... 

4 \iT57j 6rgS 



(3.11) 



At 12 days the observed OSSE flux at 100 keV (Leising et al. 1994) corresponded 
to 4.3 X 10 3 ^ erg s - ^ keV - ^. With V4. = 2 the required mass loss rates are 
M_$/v w i = 14.6, 9.8, and 8.2 for s = 2 and T e g = 1, 2, and 3, respectively. 
With s = 1.5 we get M_^/v w i = 12 A, 8.5, and 7.1, respectively, at 10^ cm. These 
are considerably higher than those derived above from the ROSAT and ASCA fluxes 
using equation (3.10). 

The ratio of the OSSE and ROSAT luminosities is given by equations (3.10) 
and (3.11), 

T cs (100 keV) ( 3 -3)(4-3) 2 T e r 24 , 8 (1 + 3 T e cs9 ) _ L16/Te cs9 

L rev (lkeV) • 4(n-3)(n-4)2 T° c 5 s9 

(3-12) 

The observed ratio is T(100 keV) /L(l keV) ~ 0.33. The theoretical ratio is sensitive 
to both n and T e cs . However, for s = 2, n ^ 6 and 2 ^ T e rev $ ^ 5 the maximum 
ratio is only 0.11 for T e cs ^ 3 X 10^ K. With s = 1.5 the maximum ratio is 0.25. 
The strong dependence on n comes from the ratio of the emission measures of the 
reverse and circumstellar shocks, //both shocks are adiabatic the observations thus 
favor low values of both n and s. 

To include cooling and lack of ion-electron equipartition, as well as to include 
a non-uniform density, we have calculated the shock structure numerically. In the 
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adiabatic case we have varied the power law index from n = 5 to n = 7. The 
mass loss rate has been adjusted to give approximate agreement with the ROSAT 
and ASCA luminosities at ~ 10 days. Except where otherwise stated, the ejecta 
velocity at 10 days has been set at ^(lO days) 2. The initial radius, R 0} is taken 
as the radius where the ejecta and wind densities are equal; for M_§/v w \ = 2, 
i?o~5x 10 13 cm. 

In figure 3 we show for T e = Tcoul the electron and ion temperatures and the 
density in the interaction region for n = 6 and M-^/v w \ = 3, at day 10. The 
density structure is similar to the similarity solution with a dense, shocked ejecta 
shell and a low density shocked wind. The only difference is close to the contact 
discontinuity, where the finite initial radius gives a less peaked density. Immediately 
behind the reverse shock the electron temperature is a factor ~ 2 smaller than the 
ion temperature, but has reached the ion temperature halfway through the shell. 
The main result of the numerical calculations is that behind the wind shock the 
electron temperature is much lower than the ion temperature, and T e does not rise 
much higher than ~ 10 9 K at any time. This has important consequences for the 
X-ray emission from the wind shock. 

For M_§/v w i = 3 we find cooling to become important for n ^ 8. This leads to 
a catastrophic fall in T e below ~ 2 X 10^ K, down to ~ 5 X 10 3 K where photoelec- 
tric and Compton heating stabilizes the temperature (CF94). These calculations 
basically confirm the estimates from equation (3.7), although the exact borderline 
between the adiabatic and cooling case is constrained to n ^ 7 for M_§/v w \ ^ 2 
and ^(lO days) = 2. A higher value of V , or lower M, relaxes this condition (eqn. 
[3.7]). 

In the adiabatic models, most of the observed X-ray emission below ~ 10 keV 
comes from the reverse shock wave, while the forward shock front is a source of 
hard X-rays ( ^ 100 keV). The larger emission measure behind the reverse shock 
causes the emission from this region to dominate up to high energies. For n < 5 
the temperature of the reverse shock is comparable to that behind the circumstellar 
shock, ~ 10 9 K, if T e = T Coul . 

In Table 2 we give the luminosities at 10 days in the ROSAT, ASCA and OSSE 
bands for n = 5 - 7, ^(lO days) 2, s = 2, and for different M_^/v w \. The 
n = 7, M_$/v w i = 1.5 model gives a fair agreement with the observed fluxes 
in the ROSAT and ASCA bands. The ASCA/ROSAT ratio is uncomfortably low 
because of the low temperature in the shocked ejecta, especially close to the contact 
discontinuity. The analytical estimates from equations (3.10 - 3.12) are accurate 
to within a factor of ~ 2, and these expressions are therefore useful for estimating 
the luminosities. The most serious objection to this model for the early X-rays is 
that the OSSE luminosity is far too low, even with collisionless heating included, as 
shown by the T e = T{ on models. A larger M increases the OSSE luminosity, but then 
overproduces the low energy luminosities. Increasing V(10 days) also increases the 
OSSE flux, but even models with ^(lO days) = 3 are not sufficient for this M and 
n. Models with n = 6, or lower, in combination with a larger M give an improved 
fit to the combined X-ray and OSSE observations. The lower n results in a higher 
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reverse shock temperature, and larger 50 - 150 keV flux, while a higher M is needed 
to compensate for the decrease in emission measure with decreasing n (eqn. [3.10]). 
Because of the higher ion temperature, electron-ion equipartition is slower, and we 
find that, unless collisionless heating is important, it is impossible to reproduce the 
OSSE luminosity for the n = 6 model. Retaining the other parameters of the model, 
equipartition increases the calculated OSSE luminosity by a factor ~ 3 compared to 
the T e = Tcoul model, still a factor ~ 4 low compared to the observations. Models 
with n = 5 in Table 2 continue this trend, but these also fail to give a high enough 
OSSE flux, confirming the analytical discussion. 

The spectrum of the model in figure 3 at 10 days is shown in figure 4, together 
with the T e = T{ on model. The contributions from the reverse and circumstellar 
shocks are shown separately, and it is clear that the circumstellar contribution is 
small even in the OSSE range in the pure Coulomb heating case. The circumstellar 
shock is important only for the T e = T{ on model, and then only above ~ 100 
keV. In these spectra we have not included the Comptonization of the photospheric 
photons. Because of the very high electron temperature in the T e = Tj 0n model, 
T e 6.9 X 10^ K, Comptonization dominates the X-ray emission below ~ 10 keV 
in contradiction to observations (§3.2 and Table 4). 

If the wind density decreases slower than r , as indicated by the radio obser- 
vations during the first ~ 100 days, this leads to a slower decline of the X-ray flux. 
If the reverse shock is adiabatic at the observed epochs it will remain so. The total 
X-ray luminosity decreases as £ - [12-7s+n(2s-3)]/(n-s)_ reverse shock temper- 

ature is falling as t~^^~ s ^^ n ~ s \ and the luminosity falling in a given band with 
E <C kT rev depends on time as 

L(E) oc t -[(2 S -3)n-5 S +6]/(n- S )_ (3^ 

For s = 2, n = 7, L(E) oc t -0 - 6 , while s = 1.5 and n = 7 give L(E) oc t - 27 . 
Neglecting the weak temperature dependence on the Gaunt factor, the observed 
decline could in the adiabatic case be reproduced by s ~ 1.85. However, we again 
stress that this case is probably not applicable to SN 1993J. 

3.2. The radiative case 

Radiative cooling at the reverse shock front becomes important for n ^ 8. 
Above 2 X 10 7 K, the cooling is dominated by bremsstrahlung and the cooling 
time is given by equation (3.7). Below 2 X 10 7 K, the cooling is dominated by lines 
and the cooling rate is A(T e ) 3.4 X 10 -23 T-T ' 67 erg cm 3 s _1 (Raymond, Cox & 
Smith 1976). In the radiative case, the cooling time of the reverse shock is therefore 



3.5 x io», < 4 : 3 " 3 - 3 ' 43 ' v?-*- (M^V 1 (-L 



s. 



cod ~ (n- 3)(n - 4)(n - s) 3 - 34 4 V "wl J \ll-57 days, 

(3-14) 

The reverse shock is radiative up to at least one year for all M_5 v ^ ^ 2, V4 ^ 2 
and n ^ 15. However, the cooling time is sensitive to s and especially to n, 
^cool ^ n - ^' 34 for n ^> 4. A shallower density gradient, as is likely in the interior, 
increases the cooling time scale considerably. 
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All X-rays from the reverse shock are absorbed before ~ 50 days, and only 
emission from the circumstellar shock is observed. The cool shell between the 
reverse shock and the contact surface has a column density 



,21 2 (n - 4) / M_ 5 \ x _ s / t_ 

(3-s)(A-s)\v wl ) 4 ^11.57 days, 



N cool = 9.07 x 10^ /o \ v<| ( ) V/"* ( r^-^ ) cm"'. (3.15) 

X-rays below 



E T= i = 1.23 



2 (n - 4) M_ 5 



(3 - s)(4 - s) v w i 



3/8 ,r 1VS / i x -3( s -l)/8 

-3(,-l)/8 1 \ ^ 

4 V 11.57 days/ 



(3.16) 

are absorbed, and the temperature of the reverse shock therefore has to be less 
than ~ E T= i/k. At 10 days this is for SN 1993J satisfied for n ^ 13 if V\ ^ 2 and 
M_$/v w i sa 4. The late X-ray observations, discussed below, indicate n > 20 and 
the cooling constraint is well satisfied. 

If cooling, the total energy emitted from the reverse shock is 
(3 - sf{n - 3)(n - 4) MV 3 ( r \ 2 -s 



2 (4 — s)(n — s) A v w Vr, 



o 



41 2 (3 - s) 2 (n - 3)(n - 4) M . 
1.57 x 10 41 - V ; V A ; 5 




(4-s)(n-s) 3 v w i 

(3-17) 

Half of this is emitted outwards. Both the energy absorbed by the ejecta and the 
cool shell contribute to the bolometric luminosity of the supernova. As long as 
cooling is important, the luminosity decreases as L rev oc £ - (15-6s+sn-2n)/(n-s) _ 
This energy input is mainly deposited in the outer ejecta close to the shock, and in 
the cool, shocked gas. This may give rise to broad emission lines in the UV and to 
Ha (e.g., Fransson 1984; CF94). A high density in the ejecta and cooling shell at 
early times may thermalize the lines (Baron et al. 1994). 

At 10 9 K < T e < 5 x 10 9 K the Gaunt factor at 1 keV is g ff = 3.45 T g - 663 to 
within ~ 2%. The ROSAT luminosity is then 

Ml ^V) « 4 ' 1 X 1037 ( Tg ' 16 e -0-0116/T 9 



(3-3) 




t 



3-2s 



erg s ^ keV ^. 



V v wl J \ 1 1.57 days J 

(3.18) 

With V4 = 2 and ( = 0.86 we obtain M_ 5 /u wl 6.7 (s = 2) and M_ 5 /u wl 7.5 
(s = 1.5), nearly independent on temperature. The ASCA/ROSAT ratio is 2.3 - 
2.6 for 10 9 K < T e < 5 X 10 9 K, within the observed range. 

Also in this case, the OSSE flux is given by equation (3.11), and the same values 
for M /v W} depending on temperature, are required. The OSSE/ROSAT ratio is just 
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given by the form of the free-free spectrum and can be approximated by 

L ( 100 keV ) _ ff//(100 keV) i.i 6/Te 9 ^ 0.143 _-1.16/T e 9 ( . m 

L(l keV) " g ff (l keV) ~ e 9 ' 

For T e ^ 1.2 x 10 9 K the ratio is within a factor two of the observed one, which 
is L(100 keV)/L(l keV) ~ 0.33. This ratio is considerably easier to obtain in the 
cooling scenario. 

The importance of cooling, as well as departure from electron-ion equipartition, 
makes a numerical calculation necessary. For the low reverse shock temperatures in 
these models the cooling time is very short, with gas cooling from ~ 10^ K down to 
~ 10^ K in less than one day at an age of 10 days for n <; 15 (eq. [3.14]). The short 
cooling time behind the reverse shock makes an accurate calculation of its structure 
difficult. To tackle this problem we have used two approaches, either by a modified 
hydrodynamic method, or by using an extension of the similarity solution. 

In the hydrodynamic approach we avoid too short time steps by merging the 
cool shells continously. Even then the structure was difficult to resolve, and the 
total number of time steps became more than 10^. Instead of attempting even 
finer zoning and smaller time steps, we take advantage of the short cooling time 
and calculate the structure of the reverse shock assuming a stationary structure, 
with realistic cooling and emission from the gas. The temperature and density 
immediately behind the shock are determined from the time dependent solution. 

In the other approach we use the similarity solution for an adiabatic circum- 
stellar shock and a radiative reverse shock (CF94). This applies if cooling of the 
circumstellar shock is negligible and if the finite initial radius is ignored, i.e. after 
t 3> Ro/Vg ~ 0.2 days. The electron-ion heating of the shocked circumstellar gas 
is calculated by following each mass element, starting at the shock. The structure 
of the cooling shock is calculated as above. After solving the ordinary differential 
equations describing the density and temperature behind the reverse shock and in- 
cluding cooling, we calculate the spectrum from the cooling gas. The emission is 
assumed to be in equilibrium at the local temperature and density. The emissivity 
is calculated using recombination and collisional ionization rates from Arnaud & 
Rothenflug (1985), and collisional excitation rates from Gaetz & Salpeter (1983), 
Mewe & Gronenschild (1981), and Raymond & Smith (1979). To calculate the pho- 
toelectric absorption of the cool gas between the shocks we assume this gas to be 
neutral. The spectrum from the circumstellar shock is treated as for the adiabatic 
shock. 

Figure 5 shows the temperature and density at 10 days for the s = 2, n = 25, 
M_$/v w i = 3.5 and T e = T(j ou \ model, calculated with the full hydrodynamics. The 
mass loss rate was adjusted to approximately give the observed ROSAT and ASCA 
luminosities at 10 days. Both density and temperature are nearly constant behind 
the circumstellar shock. As with the model in figure 3, the electron temperature 
is only ~ 1.6 X 10 9 K, much lower than the ion temperature. The cooling shell 
behind the reverse shock is very thin and can only be seen as a spike in the density 
distribution and a similar dip in the temperature. In figure 6 we show an expanded 
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view of the structure of the reverse shock, calculated from the stationary solution. 
The cooling takes place over a very narrow spatial interval, and has basically the 
structure of a thermal instability. Most of the energy is released in this region. The 
cool shell is likely to be broken up by instabilities in a multi-dimensional calculation. 
The column density of this shell corresponds well to that obtained from equation 
(3.15). The T e = T{ on model is similar, but has an electron temperature ~ 5 X 10^ 
K. 

In Table 3 the ROSAT and ASCA luminosities are given, as well as the OSSE 
100 keV luminosity, at 10 days for different values of s, M and assumptions about 
electron-ion equipartition. The value of n is of minor importance at this early 
epoch, as no emission from the reverse shock penetrates the cool shell. As expected 
from equation (3.18) the luminosity scales roug hly as M 2 . The mass loss rate was 
adjusted to give approximately the correct ROSAT and ASCA luminosities, and 
does not vary much with various assumptions of equipartition, n or s. We find that 
M_$/v w i (4 — 5) is required. Compared to the values of M/v w derived from 
equation (3.18), this is smaller by a factor ~ 1.7, which is mainly a result of the 
nonuniform density structure. As shown by the T e = T{ on and T e = T(j ou \ models 
the OSSE flux is sensitive to the degree of electron-ion equipartition. The T e = Tj 0n 
model has a flux higher by a factor of ~ 2, compared to the T e = T(j ou \ model. The 
OSSE/ROSAT ratio is close to the maximum value inferred from equation (3.19). 
Although somewhat low, even the T e = Tcoul models give acceptable fits to the 
observed L(100 keV) 4.3 x 10 38 erg s _1 keV" 1 . 

A high T e , as in the T e = T{ on models with collisionless heating, increases the 
Comptonized flux in the keV range considerably. The electron scattering optical 
depth behind the circumstellar shock is 



2.0 x 10 -2 > M_k i s / t , , 



In each scattering the photons increase their energy by a factor Av /v ~ 4kT e / m e c 

1. The probability of being scattered N times is . The multiple scattering 
creates a power law continuum that may reach as far up in energy as the X-ray 
regime, the slope mainly being determined by r e and T e (Fransson 1982, 1984). We 
calculate this component by a Monte Carlo method, which takes relativistic effects, 
spherical geometry and the actual shock structure into account (cf. Lundqvist & 
Fransson 1988). On day 7, r e = 0.041 - 0.056 for M_ b /v wi = 4, n(He)/n(H) = 0.1 
and 1.5 ^ s ^ 2. In Table 4 we give the slope of the spectrum and the Compton to 
free-free ratio of the luminosity in the ROSAT range for r e = 0.05, s = 2 and for 
different values of T e . To estimate the free-free luminosity we use equation (3.18) 
with n(He)/n(H) = 0.1. For low values of T e most of the Comptonized flux falls 
below 300 eV where the interstellar optical depth is greater than unity, assuming a 
galactic column density of ~ 4.3 X 10 2 ^ cm -2 (Zimmermann et al. 1994). Therefore, 
we have listed in Table 4 the Compton to free-free ratio corrected for this column 
density. For other values of s the Compton to free-free ratio scales as oc (3 — 5) - ^. 
Regardless of this, it is clear that for electron temperatures higher than ~ 1.5 X 10^ 
K, the Comptonized flux dominates the observable flux in the ROSAT range. This 
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is clearly incompatible with the flat spectrum observed by ROSAT (Zimmermann 
et al. 1994). In figure 7 we show the full spectrum, including Comptonization, 
for T e = T ion w 5 x 10 9 K and T e = T Cou l ~ 1-5 x 10 9 K on day 10, clearly 
demonstrating the difference. From the absence of X-rays due to Comptonization 
one can therefore rule out equipartition in the circumstellar shock. 

The luminosity from the reverse shock is high. For a mass loss rate of 
M_$/v w i = 4 we get luminosities in the range (1.9 - 3.6) x 10 41 erg s" 1 at 10 
days, for s = 1.5 — 2.0 and n 25 — 30. On day 12 the bolometric luminosity was 
~ 2 x 10 42 erg s _1 (Lewis et al. 1994; Richmond et al. 1994). 

The most characteristic aspect of the radiative scenario is the evolution of 
the X-ray flux. Because electron-ion equipartition may not be achieved for the 
circumstellar shock, and because relativistic effects are important for the free-free 
emission, the time evolution of the X-rays differs from that given by equation (3.13). 
As long as the cool shell is optically thick the spectral luminosity evolves as given 
by equation (3.18), 

L(E) oc T e - 16 *(3-2s)(n-3)/(n-s)_ 

For n ^> 10, L(E) oc t~^ s ~^\ so that for a circumstellar density gradient approach- 
ing s = 1.5 the luminosities in the ROSAT and ASCA bands are nearly constant 
with time. The observed slow decline up to ~ 40 days, L(E) oc t _ ^-30 (Zimmermann 
et al. 1994), is reproduced with s ~ 1.65. 

After ~ 100 days the soft X-ray flux from the reverse shock begins to penetrate 
the cool shell. The column density of the absorbing shell decreases as N coo i oc 
(n— 3)(s— l)/(n— s) ^ ag } on g as the reV erse shock is radiative. For n ^> 10, N coo i oc 
H* -1 ). As t(E w kTrev) b ecomes less than unity the total soft X-ray luminosity 
evolves according to equation (3.17). When t ^ t coo i, with t coo \ being given by 
equation (3.7), the flux decreases. 

An important piece of information in favor of the cooling scenario came from 
the ROSAT observation from November 1993 (Zimmermann et al. 1993; 1994). 
The temperature of the X-ray emission on November 1-2 1993 was found to be 
only ~ 0.5 keV. This estimate is uncertain because line emission dominates below 
~ 2 X 10^ K. Here we propose that this emission originates from the reverse shock, 
requiring n ^ [15.1 (3 — s) V4 + s]. This implies a very steep gradient; for V4 ~ 1.5, 
which is reasonable at 220 days, we get n 24.6 for s = 2 and n 32.2 for 
s = 1.65. Equation (3.14) shows that the reverse shock is radiative at 220 days for 
all M_ 5 v~\ > 2, V4 < 2 and n > 15. 

In November 1993, the expected column density of the cool gas behind the shock 

is 

N cool = 2.6 x 10 20 {U ~ 4) M - 5 cm" 2 , (3.22) 

for s = 2, and a factor 2.7 higher for s = 1.5. In both cases V4 = 1.5 has been 
assumed. Zimmermann et al. (1994) find N coo i fti 6 X 10 2 ^ cm -2 . For n ^ 15 
and M_5 ^ 2 equation (3.22) gives N coo i ^ 1 X 10 22 cm -2 . The s = 1.5 case 
is a factor ~ 3 higher. Both estimates are considerably higher than is observed. 
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However, Chevalier, Blondin & Emmering (1992) find that the interaction region is 
unstable, and it is likely that the absorption column density can be appreciably lower 
than estimated from equation (3.22). Typically, Chevalier et al. find a factor ~ 3 
variation in the column density for different lines of sight. Also, the circumstellar 
density may decrease faster at large radii, compared to the density law (see 

§2). 

On day 220 the total luminosity from the reverse shock is approximately given 

by 

f 7.3 x 10 39 V? erg s" 1 (s = 2) / N 

5.0 x 10 40 ^ V?- b erg s" 1 (s = 1.5), 

where we have taken n = 20 in both cases. L rev is fairly insensitive to n, but 
sensitive to s. We estimate the observed luminosity by Zimmermann et al. to be 
~ 1 X 10 39 erg s~\ lower than from equation (3.23) for reasonable M_5 and V . 
Because of the high column density a large fraction of the luminosity is absorbed 
(see below), so the discrepancy is expected. Much of the remaining luminosity 
may be emitted as high velocity optical and UV lines from the ejecta and the cool, 
absorbing shell. Especially in the s = 1.5 case the soft X-ray luminosity may be 
disturbingly high compared to the bolometric luminosity. 

Figure 8 shows the evolution of the calculated ROSAT, ASCA and OSSE 
luminosities for the s = 2, n = 25, M-^/v w \ = 3.5, T e = Tcoul case. While 
the OSSE flux decreases steadily with time, especially the ROSAT flux undergoes 
a dramatic change at ~ 65 days, shown in figure 8 as a minimum and then a slight 
increase of the light curve. The flattening is a result of the gradually increasing 
transparency of the cool shell: Njj oc t~^ for n ^> 10. The flux from the reverse 
shock, as shown by the dashed line, contributes an increasing fraction to the total 
flux in the ROSAT band. 

The spectra of the T e = T(j ou \ and T e = T{ on models at 10 days and 200 
days are shown in figure 9. While the 10 day spectrum is a smooth, free-free 
one-temperature spectrum from the circumstellar shock, the 200 day spectrum is 
completely dominated by the soft X-rays from the reverse shock. A large number of 
lines are seen, the strongest due to Fe XI-XIII and Ne X at ~ 0.91 — 1.02 keV, Mg 
XI-XII at ~ 1.35-1.48 keV and Si XIII at ~ 1.87 keV. Below ~ 1 keV photoelectric 
absorption due to the oxygen K-shell and iron L-shell gives an abrupt cutoff. The 
cutoff energy decreases with time according to equation (3.16). After the cool shell 
has become transparent the calculated ROSAT flux above the cutoff energy remains 
approximately constant, as long as the reverse shock is radiative and the density 
gradient is constant. We also note the increased extension of the high energy tail 
at 10 days for the T e = T{ on case. 

As with the adiabatic models, the calculated ROSAT luminosity declines too 
fast for the s = 2 model, and a less steep density gradient is indicated. Figure 10 
shows the light curves for p oc r - ^, n = 30 and T e = T(j ou }, and an improved 
agreement with observations over the first 100 days is found. Because of the 
higher temperature in the s = 1.7 model the flux in the OSSE range has increased 
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considerably, compared to that in the s = 2 model. 

3.3. Summary of hydrodynamics 

We favor the radiative case, compared to the adiabatic case, for several reasons. 
The density gradient of the ejecta does not need to be uncomfortably small, as is 
the case for the adiabatic scenario. The value of n needed in the radiative case is 
close to that required to reproduce the radio light curves (§ 2), and the lower limit 
to n implied by the VLBI observations (Bartel et al. 1994). The ROSAT/OSSE 
ratio is considerably easier to reproduce, and the mass loss rate derived from the 
X-rays is approximately the same as that derived from radio. In the adiabatic case 
the discrepancy is at least a factor two. To reproduce the OSSE flux the adiabatic 
models require substantial collisionless heating. The absence of Comptonized X- 
rays, however, rules out equipartition in the circumstellar shock, because this would 
give T e <; 5x10^ K. Finally, it is extremely difficult to explain the soft X-rays at 220 
days in the adiabatic model, while they are a natural consequence of the radiative 
model. Possible problems with the cooling scenario may be the high bolometric 
luminosity, and the absence of broad emission lines. However, as noted in §2, there 
may be a transition to a steeper density gradient in the circumstellar medium at 
~ (2 — 3) X 10^ cm. This would reduce the luminosity from the reverse shock (eqn. 
[3.23]) after 100 - 150 days, thus making it more likely to be compatible to the 
observed bolometric luminosity. 

Similar calculations for the X-ray emission to ours have been done by Suzuki 
et al. (1993). Their conclusions that a high velocity and shallow gradient are 
necessary differ from ours. The main reason for this is that they only consider 
the X-ray observations during the first months. From our discussion above, the 
adiabatic and radiative cases do not differ substantially at these epochs, although 
the OSSE/ASCA ratio is more easily satisfied in the radiative case. However, a 
shallow density gradient is difficult to invoke in order to explain the November X- 
ray observation, and would probably also have problems with avoiding too high an 
X-ray flux due to Comptonization during the first weeks. The mass density Suzuki 
et al. use, M_^/v w \ ^ 1, is a factor <; 4 lower than we deduce from the radio 
observations, whereas the initial shock radius they use, 2 X 10^ cm, is considerably 
larger than ours, ~ 5 X 10^ cm. Both a higher mass loss rate and a smaller initial 
radius make cooling important, and we find that the n = 12 model by Suzuki et 
al. is incompatible with the absence of X-ray absorption (and marginally also the 
n = 8 model). The general features of the light curve, as well as the relative ratio 
between the calculated ROSAT and ASCA fluxes, agree. 



4. TEMPERATURE AND IONIZATION OF THE CIRCUMSTELLAR MEDIUM 

The temperature of the circumstellar matter is crucial for estimating the mass 
loss rate from the radio observations, and the ionization structure must be known 
to calculate the absorption of the X-rays. We are also interested in the formation of 
emission and absorption lines in the gas. To calculate this we use the same method 
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as was used earlier to model the circumstellar structure of SN 1979C and SN 1980K, 
as well as the ring structure around SN 1987A (Lundqvist & Fransson 1988, 1991). 

4.1. Evolution of temperature and ionization 

Ionizing radiation is emitted during breakout of the supernova shock, as well 
as from the gas enclosed by the reverse and circumstellar shocks. Calculations 
of the radiation accompanying the shock breakout for red supergiants with radii 
(3 - 8) X 10 13 cm have been done by Falk (1978) and Klein k Chevalier (1978), who 
found an effective temperature immediately after the breakout of (1 — 2) X 10^ K, 
with a total energy of ~ 3 X 10^ ergs in a few hours. Specific calculations for SN 
1993J, assuming LTE and no scattering effects, by Shigeyama et al. (1994) result 
in a peak temperature ~ 2.9 X 10^ K and a total energy of ~ 10^ ergs. After 
~ 2 days the radiation temperature has decreased to below ~ 2 X 10^ K, and the 
ionizing effects are small. Thereafter, the photospheric effective temperature of this 
model is roughly in accord with IUE observations during the first days (Fransson 
& Sonneborn 1994) and ground-based observations (e.g., Lewis et al. 1994) at later 
times. As soon as the shocks form the ionizing radiation is dominated by Compton 
upscattering of the photospheric photons by the circumstellar shock wave (see §3.2). 
The circumstellar structure is calculated using a time dependent photoionization 
code, including all interesting elements and ionization stages. Recent updates to 
this code are summarized in Lundqvist, Fransson & Blondin (1994). Initially, the 
optical depth through the circumstellar gas at 13.6 eV is ~ 10^, too large to be 
handled by the full code. During this period, recombination can be ignored, and 
we follow the ionization using an adaptive mesh method similar to that described 
in Lundqvist & Fransson (1991). 

The recombination time of hydrogen in the circumstellar gas (assuming T = 10^ 
K) is t rec ~ 3.5 MZ5 v w \ r|g days, where r\^ is the radius in 10^ cm. If ionized only 

by the radiation from the supernova breakout, the gas inside a radius of ~ 10^ cm 
has time to recombine before the shock hits the gas. In this case, line emission from 
ions like N V, Fe X and Fe XIV is expected to be present during this period. The 
observed light curves of the lines are smeared out by the finite light traveling time, 
t ec f l0 = 2 r/c 0.77 r\^ days. However, the circumstellar gas is prevented from 
recombining by the radiation from the interaction region if the circumstellar and 
reverse shocks form prior to t rec . Exactly when shock formation occurs is uncertain 
(Epstein 1981). If preacceleration of the gas by the radiation accompanying the 
shock breakout is efficient (see §4.4), a viscous shock may only develop after about 
one expansion time. Otherwise, it may form earlier. The optical depth to electron 
scattering of the circumstellar gas from radius R out to infinity for s = 2 is unity 
for R Te =l ~ 2 X 10 M_5 v w ^ cm. The scattering photosphere may therefore be 
located at a considerable distance out in the circumstellar gas. This delays both 
the escape of the radiation and the formation of the viscous shock. In the case of 
s = 1.5, or a cavity around the progenitor formed by the influence of a binary, the 
delay of the burst is not important. 

We have calculated the ionization and temperature structure of the circumstellar 
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gas for two models: one has M-^/v w \ = 5, s = 1.7 and n = 30, and the other 
M_$/v w i = 4, s = 2 and n = 25. Both give good fits to the X-ray emission as 
discussed in §3.2. We adopt n(He)/n(H) = 0.3, n(N)/n(0) = 1.5 and n(N)/n(C) = 
3.5 as was found for the circumstellar gas around SN 1987A (Lundqvist et al. 1994). 
Other elements have solar abundances, and the overall metallicity is also solar. For 
the first two days we take the evolution of the photosphere from Shigeyama et 
al. (1994). Because the color temperature is likely to be larger than the effective 
temperature, we have multiplied the latter by a correction factor q = T C0 \ m .jT e j j , 
retaining the bolometric luminosity in the model by Shigeyama et al. . Klein & 
Chevalier (1978) found in their models that q ^ 2, and Ensman & Burrows (1992) 
found q ~ (1.0 — 1.8) for SN 1987A. Here we use q = 1.35. We also include the time 
delay between the ionizing photons from the limb of the supernova and its center. 
For times later than two days we have joined the model by Shigeyama et al. to IUE 
(Fransson & Sonneborn 1994) and ground-based observations (Lewis et al. 1994). 
The radiation from the circumstellar and reverse shocks is included as discussed 
above, assuming the electrons of the gas behind the circumstellar shock are heated 
by Coulomb collisions only. Furthermore, we assume that the viscous shock forms 
within the first expansion time scale, i.e., prior to ~ 0.23 days. The exact time of 
shock formation is not important for the temperature and ionization structures of 
the circumstellar gas at observable times as long as the viscous shock does not form 
later than after ~ 0.5 days (see below). 

During the supernova breakout, hydrogen and helium in the circumstellar gas 
quickly become fully ionized; even at 10^ cm the ionization takes less than 100 
s. Carbon, nitrogen and oxygen all attain their He-like stages, and iron close to 
the supernova is ionized up to Fe XVII. The temperature of the circumstellar gas, 
T^windi becomes fairly constant, ~ 3 X 10^ K, throughout the wind, scaling roughly 
linearly with the adopted value of q. T w i n d is also a weak function of the He/H- 
ratio; for n(He)/n(H) = 0.1 instead of 0.3, T w ^ nc i is reduced to ~ 2.8 X 10^ K. As 
the color temperature of the burst decreases, the gas close to the supernova cools by 
Compton cooling. Further out in the circumstellar gas the initial structure is more 
or less frozen-in. If the viscous shock forms after ~ 0.5 days, the photospheric flux 
above 100 eV is low enough for recombination to N V close to the ejecta to become 
important. In our models, shock formation occurs earlier and recombination is 
negligible; the circumstellar gas becomes completely ionized by the radiation from 
the shocks, once these have formed. At the same time, T w i n( i increases rapidly. 

In figure 11 we show the evolution of T w i n d at 1.5 X R S} which is roughly the 
temperature to be used in equation (2.3). The peak is due to Compton heating of 
the circumstellar gas by the radiation from the circumstellar shock. This radiation is 
during the first days totally dominated by Comptonization of photospheric photons, 
and completely ionizes the circumstellar gas. The flux due to Comptonization 
decreases rapidly with time, and as a result so does also the temperature of the 
circumstellar gas close to the shock. After ~ 10 — 15 days, the temperature remains 
roughly constant at a level slightly in excess of 10^ K. The ionization structure 
evolves less dramatically; around day 5, the gas is still essentially fully ionized, 
except for iron which is dominated by Fe XIV. After ~ 5 days, high populations 
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of hydrogen-like C, N and start to build up, especially in the s = 1.7 model. 
This has important effects on the circumstellar X-ray absorption (see §4.2). After 
~ 30 days the degree of ionization stays roughly constant. At 200 days, the s = 1.7 
model is dominated by roughly equal populations of fully ionized and hydrogen-like 
stages of carbon and nitrogen, whereas oxygen is mainly in O VII-VIII and iron in 
Fe VIII-IX. The gas in the s = 2 model has recombined even less. An important 
difference between the two models is that the zone with partial recombination is 
more extended in the s = 1.7 model, reaching out to ~ 15 X R S} whereas in the 
s = 2 model it is mainly gas inside 5 X R s that is affected by recombination. 

4.2. Circumstellar X-ray absorption 

The column density of the circumstellar gas is for n(He)/n(H) = 0.1 and £ = 1.3 
(§2) given by 

N(H) CS = 2.1 x 10 22 , M ~ 5 V}- 8 ( *—) cm" 2 . (4.1) 

V ; (s-l)v wl 4 \8. 90 days/ V ; 

ROSAT observations show the total column density to be consistent with the 
galactic absorption in the M 81 direction, ~ 4.3 X 10 2 ^ cm -2 (Zimmermann et 
al. 1994). With t w 10 days and V 4 « 2 we find M_ b /v wi < 0.045 (0.025) for 
s = 2.0 (1.7). These values are clearly incompatible with the value of M/v w 
derived from the radio and X-ray observations, unless the circumstellar gas is almost 
completely ionized so that K-shell electrons of C, N and O, and L-shell electrons of 
Fe are removed. In figure 12 we show the optical depth through the circumstellar 
gas at 1 keV, r(l keV), as a function of time for the s = 1.7 and s = 2.0 models. It is 
seen that with the fading of the radiation due to Comptonization, r(l keV) increases 
rapidly and peaks around day 10 — 20 at r(l keV) ~ 0.3 (1.5) for s = 2.0 (1.7). 
This is higher than r(l keV) ~ 0.1 inferred from N(H) ~ 4.3 X 10 2 ^ cm -2 . The 
high value we find for the s = 1.7 model can be suppressed if the overall metallicity 
is lower than solar, the n(N)/n(0) value is higher than 1.5, or if the p w oc r~^^ 
dependence steepens at large radii, as may be indicated by observations at 3 mm 
(Phillips & Kulkarni 1993b). r(l keV) may also decrease if the circumstellar gas is 
non-spherically symmetric. 

As was discussed in §4.1, the temperature and ionization structures stay roughly 
constant after ~ 20 — 30 days. As a result, the X-ray optical depth through the 
circumstellar gas decreases with time after ~ 20 days (figure 12). We therefore think 
it is unlikely that the observed increase in column density of absorbing gas up to 
~ 200 days is due to recombination of the circumstellar gas, as has been advocated 
(Zimmermann et al. 1994). As discussed in §3.2, we believe the softening of the 
X-ray spectrum is instead a result of the cool shell between the circumstellar and 
reverse shocks becoming transparent. 

4.3. Implications for the radio absorption and the optical and UV line emission 

Our calculations give very high circumstellar temperatures close to the shock 
for t ^ 10 days (fig. 11). This reduces tjj (eqn. [2.3]) so that tjj becomes less 
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than unity at 1.3 and 2 cm at all times, in contradiction to the observations. Our 
calculations also show too low a population of N V to be compatible with the IUE 
observations. These problems may be solved in either of two ways: a low-density 
region close to the supernova (cavity), decreasing the early circumstellar emission 
from the shocks, or alternatively, clumping of the circumstellar medium. 

Clumping will lead to a higher NV A 1240 luminosity than in our model, 
though only if the ionizing luminosity does not increase appreciably with increasing 
circumstellar density. If the clumps have a filling factor, / c /, and the the ionizing 
luminosity is kept at a constant level, the line emission increases roughly as oc 
f~l^ compared to that for a smooth wind, while the effect of clumping on the 
ionizing luminosity is more uncertain. If the clumps have a spherically symmetric 
distribution, and they dissolve quickly behind the circumstellar shock, the ionizing 
flux should be similar to that for a smooth wind. If the clumps do not dissolve 
efficiently, the free-free to Compton ratio should be higher. The maximum ratio 
could be <; fcl^i an d occurs for very large clumps, or a disk. From numerical 
calculations we find that when the contribution due to Comptonization is kept at 
the same level as in §4.1 a high enough N V A 1240 luminosity is produced if /„; is 
of the order 10 . 

Clumping also affects the circumstellar absorption. In particular, tjj oc p , 

which means that as long as the wind is optically thick, tjj oc ■ If the ionizing 
flux is kept fixed, tjj increases even faster with decreasing / c /, since a higher density 
results in faster cooling of the circumstellar gas, and hence an increase of the free- 
free absorption (eqn. [2.3]). The clumpy scenario may thus give rise to both N 
V zones in the circumstellar gas, and also explain the absence of detected radio 
emission during the first ~ 5 days. If f c \ is independent of radius, we still need an 
overall distribution r _s , with s < 2, to explain the shape of the radio and X-ray 
light curves. 

A clue to whether clumping was important for SN 1993J may come from the 
observed Ha flux. Using the results by Cumming et al. (1994), and E^_y = 0.18 
(Fransson & Sonneborn 1994), one finds that the dereddened flux on day 2.2 

on 1 

corresponds to an Ha luminosity of £jj a 1.2 X 10 erg s . For a spherically 
symmetric wind, this means that if the gas emitting the line resides in a region 
between R s and r, the mass loss needed to explain the Ha emission corresponds to 



m 33.1 & (2s - 3)» *T» « ( * f" (l - (^ 3 - 2A "° 5 

v wi Jcl K ' 5 V10 15 cm7 I \R 8 , 



assuming Case B emissivity (e.g., Osterbrock 1989), n(He)/n(H) = 0.3 and s > 1.5. 
The maximum radius, r, is approximately given by the light echo radius, r = ct/2. 
For T = 10 5 K, f cl = 1 and t = 2.2 days, equation (4.2) gives M_ 5 /v wi 22 - 24 
for 1.6 ^ s ^ 2. This is a factor of ~ 10 higher than used in §4.1. A similar value is 
found for N V A 1240, if the ionizing luminosity is kept at the same level as in §4.1. 

Comptonization is also suppressed if the progenitor was surrounded by a low- 
density region, since r e oc M cs = M(R S — R ) / v w (s = 2), where M cs is the swept up 
circumstellar mass and R is the radius of the low-density region. Slower Coulomb 
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heating leads to a lower electron temperature. If the low-density region extends to 
~ (0.5 — 1) X 10^ cm, the early high level of Comptonization is avoided, giving a 
lower circumstellar temperature outside this radius than in §4.1, and hence a higher 
radio absorption. In the case of a low-density region, the radio flux is lower also 

because the radio emissivity scales as oc Pw~^^^ ■ A low-density region therefore 
can explain the observed onset of the radio emission after ~ 5 days. Unless the wind 
exterior to the low-density region is clumpy, it is necessary in this scenario to invoke 
radiative excitation to explain the observed strong NV A 1240 and Ha flux. While 
the N V line should be dominated by direct radiative excitation, Ha is probably 
boosted also by excitation to the n = 3 level by L/3. We find that a very steep 
spectrum due to Comptonization is necessary for radiative excitation to operate; 
at 2.2 days, a ~ 4 — 5. This implies that if radiative excitation is important, a 
very low density close to the progenitor is necessary; the low-density region would 
essentially resemble a cavity. 

4.4. Preacceleration 

The strong burst of radiation during the first few days may have important 
dynamical effects on the circumstellar gas. Assuming that electron scattering 
dominates the opacity, an initial UV-burst of 10 4 ^-E49 erg s - ^ preaccelerates the 
circumstellar gas to 

vcs = 1.42 x 10 4 E i9 V 4 ~ 2 (^) tf ayg km s" 1 (4.3) 

where R s is the shock radius and r the distance from the supernova. With V\ ~ 2 
the preshock velocity on day 2 is therefore ~ 900 -E49 km s~\ decreasing with 
time as t~ . Observations of Ha emission on March 29, 1993 show evidence for 
velocities of ~ 1000 km s _1 (Cumming et al. 1994). Also the NV A A 1238 - 1242 
lines indicate similar velocities (Fransson & Sonneborn 1994). We have integrated 
the energy from the light curve model for SN 1993J by Shigeyama et al. (1994), 
and find up to day 2.0, E = 6 X 10 4 ^ ergs. In addition to the radiation from the 
photosphere, the circumstellar gas is also accelerated by the radiation from the 
circumstellar shock. In the s = 2.0 (1.7) model this contribution to E amounts to 
~ 5 X 10 4 ^ (~ 2 X 10 4 ^) ergs up to day 2.0. Together, the flux from the photosphere 
and that from the circumstellar shock may be sufficient to explain the observed 
preacceleration. In the adiabatic model, n is small and V oc t~^~ s ^ n ~ s ^ is large 
at very early times (eqn. [2.2]). Therefore, it is more difficult to obtain enough 
preacceleration in the adiabatic model, which provides an additional, indirect piece 
of evidence in favor of the radiative models. 

Line acceleration may also be important (Fransson 1986). Although the ef- 
ficiency requires a detailed study, it was argued that the lines can increase the 
acceleration by a factor 4 - 10, compared to electron scattering. Line acceleration 
requires that the wind is not completely ionized, which may be fullfilled initially 
inside ~ 10^ cm, and certainly has to be the case in the gas emitting the NV 
lines. Because preacceleration is strongly dependent on radius, the line widths are 
expected to decrease with time. Fransson (1986) argued that the velocity should 
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depend on radius as r~ , and the line widths therefore decrease roughly as t~ . 
If the circumstellar gas is accelerated over an extended period of time, the line 
widths may decrease somewhat more slowly. 



5. DISCUSSION 

The r _ l'^ model shown in figures 1 and 2, and the models for the X-ray emission 
discussed in §3, indicate that there is evidence for a shallower density dependence 
in the circumstellar gas than an r law. The drop in the luminosity at 3 mm after 
~ 100 — 150 days (Phillips & Kulkarni 1993b; §2) indicates that the departure from 
an dependence only refers to the region traversed by the shock between ~ 8 
and ~ 150 days. For an dependence this translates into a variation in the 

mass loss rate by a factor of ~ 4 over this region. Variations in the density on even 
shorter length scales are indicated by bumps in the optically thin part of the 2 cm 
light curve after ~ 50 and ~ 100 days (Pooley & Green 1993; Weiler et al. 1994), 
and of the 3.6 cm light curve after ~ 100 days (Weiler et al. 1994). 

A density dependence different from p w oc r could be the result of either a 
variation of the mass loss rate on a time scale ~ R s /v w = V s t/v w 700 years 
before the explosion (equivalent to a length scale of ~ 2 X 10^ cm, see §2) where R s 
is the shock radius, or a similar variation in the wind velocity. A possible origin of 
this may be connected to instabilities in the nuclear burning during or after carbon 
burning (e.g., Woosley, Weaver & Taam 1980). The conditions for this to occur 
are, however, restricted to He-cores between 2.5 — 2.8 Mq. Another possibility is 
that the mass loss is due to pulsational instabilities in the envelope. However, the 
time scale for the fundamental pulsational mode in red supergiants is typically ~ 10 
years (Fox & Wood 1982), corresponding to radio variations on a time scale of ~ 1.5 
days. Because variations on time scales less than R s /c = Vt/c t/15 are smoothed 
out by light travel delays, the shortest mass loss variations that can be detected 
at ~ 50 days therefore correspond to ~ 20 years. Alpha Orionis shows evidence 
for periodicity in its mass loss on a time scale of 5.8 years and for structure in its 
circumstellar envelope on a scale of ~ 10^ cm (Goldberg 1979). 

If the progenitor was a member of a binary system, this may modify the 
circumstellar environment considerably. Woosley et al. (1994) estimate a binary 
separation a 1.5 X 10^ cm for SN 1993J. With a total mass of ~ 20 Mq this 
corresponds to a period of about 7.1 years, and an orbital velocity of ~ 45 km s — 
If the companion star also possesses a wind this may further complicate the density 
structure. An early type companion, as may be indicated by the observations by 
Aldering et al. (1994), would have a fast wind with velocity ~ 1000 km s - ^ and 
mass loss 10 - ^ — 10 Mq yr - ^. For two early-type components Stevens, Blondin & 
Pollock (1992) find that a dense interface forms due to radiative cooling of the gas 
shocked by the colliding winds. The dense, shocked gas may be helpful in explaining 
the NV A 1240 and Ha emission discussed in §4.3. A structure similar to that in 
the calculation by Stevens et al. is also seen in a model for the interacting winds 
in symbiotics (Nussbaumer & Walder 1993), which invoke one hot star and one red 
giant, and may even better represent the situation for SN 1993J. The momentum of 
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the fast wind from the companion may dominate the slow wind from the supernova 
progenitor. For M w ~ 4 X 10 - ^ Mq yr - ^ and v w ~ 10 km s~\ this occurs for a 
mass loss rate of the companion, M ^ 4 X 10 — ^M© yr~\ assuming a wind speed 
of ~ 1000 km s~l. It is likely that the wind structure was modified also by the 
orbital motion of the system; the wind density should be higher in the orbital plane 
than orthogonal to this, producing a disk-like structure with thickness depending 
on details of the wind interaction and the orbital motions of the two stars. To 
explore fully the interaction of the ejecta with a disk-like structure, multidimensional 
hydrodynamics is necessary. However, assuming a freely expanding circumstellar 
shock, Lundqvist (1994) argues that a disk structure may explain the radio light 
curves, without invoking a time varying mass loss from the progenitor, if the disk 
thickness increases roughly as oc r^. To model accurately the right turn-ons at 
different wavelengths, and since this requires a substantial line-of-sight through the 
circumstellar gas exterior to the circumstellar shock, it is unlikely that the tentative 
disk was thin. 

In §4.3, we found evidence for that the wind close to the progenitor may have 
been clumpy, or that it was surrounded by a low-density region. Because clumping, 
a cavity, and p w oc r~ s (s < 2) may be natural consequences in the binary scenario, 
we find it plausible that the progenitor of SN 1993J was a member of a binary. 
Asymmetry may also decrease somewhat the path length of the line-of-sight through 
the circumstellar gas, thus explaining the low absorbing column density inferred 
from the ROSAT observations. However, VLBI observations during the first ~ 50 
days (Bartel et al. 1994) showed no sign of asymmetry, so we cannot rule out an 
explanation based on a time varying mass loss. Future VLBI observations may 
reveal which scenario is the most likely one. If f c i <C 1 and/or the progenitor 
was surrounded by an extended cavity, our conclusion in §3.2 that the absence of 
Comptonization during the observed epochs is evidence for T e = T(j ou \ may not 
hold true. 

6. CONCLUSIONS 

Here we summarize our main conclusions: 

1. From the radio observations there are strong indications of a circumstellar 
medium with density flatter than p w oc r , at least out to ~ 2 X 10 cm. The 
observations are compatible with a radial outflow with p w oc r - ^, indicating either 
a decreasing M/v leading up to the explosion, or non-spherically symmetric mass 
loss. The mass loss rate from radio observations is in the range (2—6) X 10 — ^M© yr - ^ 
at 10^ cm for v w = 10 km s - ^. The main uncertainty is the temperature of the 
circumstellar gas. 

2. Most of the early X-ray flux is likely to originate in the shocked circumstellar 
gas. The reverse shock is radiative in our most likely scenario, and the flux from 
the postshock gas is absorbed by the cooling gas between the shock and the contact 
surface. The ROSAT and ASCA luminosities are compatible with the mass loss 
rate derived from the radio observations. The X-ray light curve indicates a flat cir- 

17.. 

cumstellar density, p w oc r , similar to that derived from the radio observations. 
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The soft X-rays observed at 220 days are explained as emission from the reverse 
shock leaking out as the shell becomes transparent. 

3. The level of the OSSE observations indicates an electron temperature of 
~ 1.5 X 10^ K. The absence of X-rays due to Comptonization gives an upper limit 
to the electron temperature of ^ 1.5 X 10^ K for the shocked circumstellar gas. 
This is a factor of ~ 4 less than the ion temperature, and indicates that collisionless 
heating in the shock front is not important. 

4. There are several indications of a steep density gradient in the ejecta. The 
turn-ons of the radio light curves, the radiative scenario for the X-rays, and the late 
soft X-rays all require n ^ 20. Additional, though weaker, evidence comes from 
preacceleration. A steep density gradient, n ^ 20, is consistent with both VLBI 
observations (Bartel et al. 1994) and with hydrodynamic models for the explosion 
(Suzuki et al. 1993). 

5. The temperature of the circumstellar medium is ~ (1 - 3) x 10 5 K for 
most of the time during the first 100 days. After the formation of the shocks 
the circumstellar medium was nearly completely ionized, in agreement with the 
observed lack of substantial X-ray absorption. To explain the long duration and 
absolute flux of the N V luminosity, as well as the flux of the Ha line, a non-uniform 
density distribution may be required. 

6. The high velocities indicated by the N V and Ha lines are explained as a 
result of preacceleration by the radiation in connection with the break out of the 
shock. Radiation from the circumstellar shock also contributes to the acceleration 
of the circumstellar gas. 

SN 1993J in many respects conforms to the same class of Type II supernovae 
as SN 1979C and SN 1980K. There are, however, indications from both radio, X- 
rays and UV that the circumstellar medium is different from a standard p w oc r 
stellar wind. Either a variable mass loss rate, or a non-radial outflow, may explain 
the density variation. The latter alternative may be expected from mass loss in a 
binary system, and is supported by the clumping and/or a cavity we infer from the 
NV A 1240 and Ha lines. 

Finally, there are several predictions for the future. Because of the expan- 
sion the column density of the cool gas absorbing the reverse X-rays decreases 
with time as in equation (3.15), and the soft X-ray cutoff is expected to de- 
crease according to equation (3.16). The temperature of the reverse shock is 
kT r ev « 1.16 (F/20000 km s" 1 ) 2 [(n - 2)/20]" 2 keV. The soft X-ray emission 
is therefore a valuable diagnostic of the density gradient and the hydrodynamic 
structure of the outer ejecta. As the column density decreases the observed X-ray 
flux from the reverse shock is expected to increase considerably, as long as the shock 
is radiative and the circumstellar density does not drop more rapidly than r . In 
the adiabatic phase it drops according to equation (3.13). 

As the radioactive input from the supernova decays, and if the circumstellar 
density does not decrease very fast, the ingoing X-ray emission is expected to 
dominate the energy input to the supernova, as is the case for such supernovae 
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as SN 1979C, SN 1980K, SN 1988Z and SN 1986J (CF94). This results in UV and 
optical emission lines mainly from the outer parts of the ejecta. The line widths 
are therefore expected to be considerably larger, (10 — 15) X 10^ km s~\ than at 
earlier epochs. If the circumstellar medium is highly clumpy also low velocity lines 
from shocked circumstellar gas may appear. In contrast to the earlier phases, high 
ionization lines like C IV A 1550 and [O III] A A 4959 — 5007 are strong, as well 
as Lya, Ha and Mg II A 2800. This is discussed in detail in CF94. The future 
evolution of SN 1993J is therefore of obvious interest for understanding both the 
ejecta and the circumstellar gas. 
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FIGURE CAPTIONS 



1. — Radio light curves at 2 cm assuming a spectral index of 7 = 0.7 together 
with observations from Pooley & Green (1993). The full line gives the "best 
fit" for the p w oc r - ^, n = 30 model, and the dashed for the p w oc r - ^, n = 30 
model. 

2. — Radio light curves at 3 mm (dashed-dotted line), 1.3 cm (short-dashed 
line), 6 cm (solid line), and 20 cm (long-dashed line) for the p w oc r - ^, 7 = 0.7, 
n = 30 model. Observations at 3 mm are from Phillips & Kulkarni (1993a) 
and those at other wavelengths from Weiler et al. (1994) 

3. — Density, electron temperature and ion temperature for the s = 2, n = 6, 
M_$/v w i = 3 model at 10 days. Only Coulomb heating of the electrons is 
assumed, leading to T e <C T{ on for the circumstellar shock. 

4. — Spectrum at 10 days for s = 2, n = 6 and M-^/v w \ = 3. The upper 
panel assumes T e = T(j ou \ and the lower T e = T{ on . The dotted lines give 
the contribution from the reverse shock, and the dashed lines that from the 
circumstellar shock. 

5. — Density, electron temperature and ion temperature for the s = 2, n = 25, 
M_$/v w i = 4.0 model at 10 days. Only Coulomb heating of the electrons is 
assumed. The density spike close to the reverse shock is due to the cool shell 
discussed in the text. 

6. — The density (dashed line) and temperature (solid line) of the radiative, 
reverse shock for the same model as in figure 5. 

7. — Spectrum at 10 days for the s = 2, n = 25, M-^/v w \ = 4.0 model 
including Comptonization of photospheric photons by the nearly relativistic 
electrons behind the circumstellar shock. Left panel is for T e = T(j ou \ and 
right for T e = T{ on . Short-dashed lines show the free-free luminosity, dotted 
lines the Comptonized luminosity and solid lines the total luminosity. Long- 
dashed lines show the observable luminosity assuming an interstellar column 
density of 4.3 x 10 /u cm _/ . With T e = T ion especially the ROSAT range is 
dominated by the Compton contribution, with a power law slope of ~ 1.0. This 
is in contradiction to the observed spectrum. 

8.— Luminosity as function of time in the ROSAT (0.1 - 2.4 keV), ASCA (1 - 
10 keV), and OSSE (50-150 keV) bands for the s = 2, n = 25, M_ 5 /v wi = 4.0, 
T e = Tcoul model. The dashed line gives the total luminosity from the reverse 
shock, while the dotted line gives the same from the circumstellar shock. The 
corresponding contributions to the ROSAT luminosity are also shown. 
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9. — Spectrum at 10 days and 200 days for the radiative reverse shock model 
with s = 2, n = 25, M-^/v w \ = 4.0. The dotted lines give the contribution 
from the reverse shock. At 10 days the hard emission from the circumstellar 
shock dominates at all energies, while at 200 days the soft component from the 
reverse shock dominates. Most of the emission from the reverse shock is due 
to line emission. Below ~ 1 keV photoelectric absorption from the cooling gas 
cuts off the reverse emission. Interstellar and circumstellar absorption are not 
included. 

10. — Same as figure 8 for s = 1.7, n = 30, M_§/v w \ = 4.5 and T e = Tfj ou i. 

11. — Temperature of the circumstellar gas at R = 1.5 X R s for s = 1.7 and 
M_$/v w i = 4.0 (dashed line), and s = 2.0 and M_§/v w \ = 5.0 (solid line). 

12. — Optical depth through the circumstellar medium at 1 keV for the two 
models in figure 11. Solid and dashed lines represent the same models as in 
figure 11. 
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TABLE 1 



Cooling times for the reverse shock 





for s = 


1.5 and 


s = 2 




n 


days 


t a 
cool 

days 


C n ,2 
days 


t a 

cool 
days 


6 


0.366 


30.39 


0.498 


11.20 


7 


0.150 


7.63 


0.199 


3.17 


8 


0.0760 


2.76 


0.0996 


1.26 


9 


0.0439 


1.25 


0.0569 


0.61 


10 


0.0277 


0.65 


0.0356 


0.34 


12 


0.0131 


0.23 


0.0166 


0.13 



t coo i is the cooling time at 1 day for M_§/v w \ = 2 and ^(lO days) = 2. 



TABLE 2 

Luminosities at 10 days in the ROSAT, ASCA and OSSE bands for s = 2. 



n 


M-5/vwl 


P0 

10" 16 g cm -3 


T e 




Te rev 

10 8 K 


Te cs 

10 8 K 


ra 

-^0. 1-2.4 


ra 

^1-10 


^100 


5 


4.0 


0.021 


Coul 


2.08 


4.5 


13.7 


2.89 


5.69 


0.61 


6 


2.0 


0.1 


Coul 


2.07 


3.1 


10.7 


2.06 


3.62 


0.16 








eq 


2.02 


3.1 


45.7 


3.07 


5.07 


0.47 


7 


1.0 


0.29 


Coul 


2.00 


2.0 


8.7 


1.77 


2.25 


0.03 




1.5 


0.44 


Coul 


2.01 


2.2 


10.5 


5.27 


5.18 


0.09 








eq 


1.95 


1.7 


51.2 


7.82 


6.70 


0.28 




2.0 


0.63 


Coul 


2.01 


2.2 


10.5 


6.50 


8.90 


0.17 



a in 10" 39 erg s 1 

h in 10 38 erg s _1 keV" 
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TABLE 3 



Luminosities at 10 days in the ROSAT, ASCA and OSSE bands 
for a radiative reverse shock. 



s 


n 


(M_ 5 /v w i) a 


T e 




Te rev 

10 8 K 


Te cs 

10 8 K 


r b 

L 0. 1-2.4 


r b 

L l-10 


^100 


1.5 


35 


4.0 


Coul 


2.00 


0.11 


18.2 


1.80 


5.58 


1.89 






5.0 


Coul 


2.00 


0.11 


19.8 


2.82 


8.72 


3.18 


1.7 


30 


4.0 


Coul 


2.00 


0.12 


17.7 


1.75 


5.43 


1.75 






5.0 


Coul 


2.00 


0.11 


19.0 


2.68 


8.28 


2.82 








eq 


2.00 


0.11 


53.7 


3.00 


9.66 


5.22 


2.0 


25 


4.0 


Coul 


2.00 


0.10 


16.6 


1.98 


5.97 


1.61 








eq 


2.00 


0.10 


58.1 


2.18 


6.72 


2.82 






5.0 


Coul 


2.00 


0.10 


18.4 


3.33 


10.0 


2.82 



a at 10 iO cm 

h in 10 39 erg s -1 

c in 10 38 erg s -1 keV~ 



TABLE 4 



Effects of Comptonization on the ROSAT emission 
at 7 days for r e = 0.05 and s = 2. 



Te cs 

10 9 K 


a 


Lcompl Lff 

N H = 


Lcompl Lff 

N H = 4.3 x 10 20 


1.0 


2.50 


1.36 


0.26 


1.5 


2.00 


8.25 


2.43 


2.0 


1.66 


25.6 


9.20 


3.0 


1.30 


101. 


49.5 


5.0 


0.96 


341. 


203. 
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Fig. 1 
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Fig. 3 

14.5 i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 



-15 



15.5 



-16 



16.5 



_-|7 I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i 

1.9xl0 15 2xl0 15 2.1xl0 15 2.2xl0 15 2.3xl0 15 2.4xl0 15 2.5xl0 15 2.6xl0 15 



Radius (cm) 



i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i 




1.9xl0 15 2xl0 15 2.1xl0 15 2.2xl0 15 2.3xl0 15 2.4xl0 15 2.5xl0 15 2.6x10 

Radius (cm) 



37 

Fig. 4 
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Fig. 5 
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Fig. 6 
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Fig. 8 
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Fig. 11 
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Fig. 12 



